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Abstract 

Spontaneous radiation in dissipative gas medium such as plasmas is investigated by 
Langevin equations and the modified Weisskopf-Wigner approximation. Since the refrac- 
tive index of gas medium is expected to be nearly unity, we shall first neglect the medium 
polarization effect. We show that absorption in plasmas may in certain case modify the 
Einstein A coefficient significantly and cause a pit in the A coefficient-density curves for 
relatively low temperature plasmas and also a pit in the A coefficient-temperature curves. 

In the next, the effect of medium polarization is taken into account in addition. To our 
surprise, its effect in certain case is quite significant. The dispersive curves show different 
behavious in different region of parameters. 


1 Introduction 

This work is motivated by an effect called by the original authors 11 ' as “quenching of Einstein A 
coefficient”: the ratio of two line intensities I(5805A) and I(312A) from a common upper level of 
Civ is reduced by an order of magnitude when the eletron density N e changes from 10 18 /cm 3 to 
10 19 /cm 3 . But there is some dispute about the interpretation of their experiments. Our work is 
to see whether this effect is possible theoretically. 

Physically, spontaneous emission is resulted from atom (ion) with “vacuum” electromagnetic 
field, therefore Einstrein A coefficient is a characteristic parameter of the total atom-field system, 
not just of atom itself. Its value may be different for different enviroment, such as cavity, or, in 
our case, the dissipative medium. 

Our approach is based on Langevin equations' 2 '. We don’t use Fermi golden rule to calculate 
the emission rate as did by Barnett et al' 3 ', because the radiating ion and radiated field do not 
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make up a closed system, their states both before and after the emission process axe not stationary 
states. 

^From Langevin equations an emission rate operator is defined, and a part of it, the sponta- 
neous emission rate operator may be separated out. 

Since the refractive index n(u) of gas medium is expected to be nearly 1, we shall first, in Sec 
II, neglect the effect of medium polarization, which will be taken into consideration in Sec III. 

2 The effect of medium absorption on Einstein A coeffi- 
cient 

Since in this section we take n(u>) = 1, both the frequencies and wave functions of light modes are 
the same as in vacuum. The reduction of Einstein A coefficient is therefore not by the alteration 
of modes of e.m. field but is caused by alteration of dynamics: the relevant operators now do 
not obey the Heisenberg equations but obey Langevin equations instead. We shall see that the 
frequency dependence of photon decay parameter will play an importance role in the present case. 

We have argued^ that the three-level ion problem may be reduced to two-level ion problem, 
the only exception is the calculation of population numbers among the levels, in which all the 
levels involved must be taken into account. 

In considering the emission of a particular ion, the whole plasmas other than the radiating ion 
will be regarded as a reservior, its effects are described by the damping and fluctuation terms in 
the Langevin equations: 

+ + ( 1 - 1 ) 

jA(t) = -^S + (i)-2ij:S3(()al tj (f)e- i l“-"">‘ + S + (t), (1.2) 

kj 

4&(f) = -r 3 [s 3 (() - S 30 ] + + h.c .] + ± 3 (t). (1.3) 

az kj 

where S± are usual atom level-changing operators, S 3 is the half of population difference operator: 
£3 = |(P 2 — A), a-kj is the photon absorption operator of mode (k ,j), j is photon polarization 
index. The free-varying phase factors have been separated out from the photon operators and 
atom level-changing oprators. 

We have defined^ an emission rate operator I(t ) as that part of — dP ^ which is caused by 
the interaction with photons, with the result 

/« = £«5k i 4i(f)S-(<)^ i( " _ " ol ‘ + he.. (2) 

kj 

Solving eq.(l.l) to get At) and substituting it into eq.(2), one may separate out the sponta- 
neous emission rate operator I sp (t ) from 7(f): 

Isp{t) = £ \g k j\ 2 f + h.c. ( 3 ) 

kj J ° 
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The Weisskopf-Wigner approximation now takes the modified formal 

£ | 5k j\ 2 e-*“~ u oUt-n-KMit-t) = + 2 i 6 u) 8 (t -t'), for t - t' > 0. (4) 

kj 

Substituting it into eq.(3), we get immediataly 

; sp = T? + ms-(o = 7ft(t). ( 5 ) 


Eq.(5), after taking the expecting value, is just the Einstein formula for spontaneous emission 
with 7 as the new A coefficient. 

Taking real part of eq.(4) and integrating it over t' from — oo to t, we get^ 


7/7 


1 

TTOJo 
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Umax 
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Uk(u) j 

Wb) 2 + K 2 (u/; 


(6) 


where 7 is the usual Einstein A coefficient in vacuum, u p is the plasma frequency and uwr denotes 
a cut-off frequency representing the limit of dipole approximation. 

We see that 7/7 is determined solely by photon decay parameter /c(a>), with no direct reference 
to atomic parameter. This is contrary to the conjectur of Aumayr et afi 5 b 

When k is small and independent of u, the right hand side of eq.(6) will reduce to 1. 
Actually, in plasmas, k is contributed by the inverse Bremstrahlung of free electrons and 
selfabsorption of ions of the same kind as the radiator. Both are u> dependent, but only the latter 
is important unless at very high plasma electron density (~ 10 21 /ott 3 )- So we shall only consider 
the latter, which will be denoted by «/( a;), in the following. 
kj(uj) is expressed by 


7rc 3 Aj 7 (Pi - P 2 ) r T 
4w§ (u - u Q ) 2 + 


(7) 


where AT/ represents the ion density, IY is the total width of the spectral line, Pi and P 2 are the 
average populations of the lower and upper levels. We have approximated the line profile by a 
Lorentzian shape in eq.(7), actually the profile is convolution of a homogeneous broadening (Stark 
broadening) and an inhomogeneous broadening (Doppler broadening). 

The frequency behavior of k/(u) is characterized by a peak at u = uq. It is this behavior 
which may cause significant reduction of 7/7. 

To explain this, we simulate the variation of kj with a; by a simple stepwise function 


ki(u) = 


Kq U ) 0 < UJ C 

k'q (jJ > Ul c 


(8) 


With Kq > Kq. 

If in the whole range ki(uj) either equals «o o r equals k' 0i the integral in the eq.(6) will amount 
to 7ru 0 . The curve a ^ ar S er height but a smaller width as compared with that of 

(oi-a/o) 2 +/c 2 ’ so S^ ves the same integral value when substituted in eq.(6). But the situation of 
stepwise function eq.(8) is different, comparing it with the above mentioned constant cases, it is 
not difficult to understand why the stepwise function will reduce the integral value. 
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One may use k(uj 0 ) ^ kj(u> 0 ) to scale the integration interval of important contribution in 
eq.(6). If kj(u>) drops down rapidly in this region, then the reduction will be large. Since Tt/2 
scales the interval of significant variation of kj(cj), it follows that tha ratio of these two scales 
IY/ 2 k/(o;o ) determines the amount of reduction of 7/7. The smaller the value of rY/ 2 /c/(u;o), the 
heavier will be the reduction. 

In this way, we see that the atomic parameters do affect the value of Einstein A coefficient, 
but in an indirect way. 

The numerical results are as follows. For line A = 5805 A, 7/7 is almost kept to be nearly 
1 in the whole range from N e = 10 15 /cro 3 to N e — 10 2O /cm 3 , because Ft/2kj(ijJq) is large, 
since the factor Pj — P2 in eq.( 7 ) is very small in this case. The result for line A = 312 A is 
shown in Fig 1 . In low density region, Tt is mainly contributed by Dopper broadening so that 
it is independent of N e , and this in turn makes k/(u>o) proportional to N e . Therefore when N e 
increases, Tt/ 2 ki(uo) becomes smaller, leading to the drop of 7/7. But when N e goes beyond a 
value about 1.5 x 10 18 /cm 3 , 7/7 turns up, because Tr is gradually dominated by Stark broadening, 
so that it is proportional to N e , which in turn leads to ,/V e -independence of kj(u 0 ). As the result, 
there is a pit or hollow in the curve 7/7 -N e . 

The dependence of 7/7 on temperature is also interesting. By a similar analysis, we have 
argued M that there will be also a pit or hollow in the curve 7/7 -T. 


3 Einstein A coefficient with index of refraction taken 
into account 


We begin our discussion by deriving the plasmas’ refractive index from microscopic equations. 
First we omit the contribution of free electrons which is less important. 

There are numerous ions in the plasmas, the effective coupling constant for ion l and e.m. field 
has a position-dependent factor e ,kx(0 . 

In gas, the orientation of atomic dipole (d)2i is random, yielding \gkj \ 2 independent of j, 




( 9 ) 


where 


u>k — kc. 


We define a collective atomic operator as usual 

WO = £e* x<,> 4>((), 

l 


( 10 ) 


the equation for dj^- ( t ) will be of the form 

j t a' kj (t) = iw t ai,(() + tG*S + , „,((). (11) 

In the case of steady-state plasmas, we may approximate the S® in the equation of <S +) kj by 
its average value |(P 2 — Pj) plus a fluctuation term. This leads to 

-f-P+kdt) = ( iu 0 - ir) 5 +k>(t) + iG k NjV(Pi - P 2 )a[At) + fluctuation term. ( 12 ) 
at 2 
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We see from eqs.(ll) and (12) that a^- and S'+.k? turn into each other. 

To derive the expression for index of refraction n, we may assume, as usually does in classical 
electrodynamics, that S +i kj will do forced vibration following d^- except a fluctuation term^l By 
this, we get an expression for 5'+,kj from eq.(12). Substituting it into eq.(ll), yields 


d. t 

dt< 


UJk + 


7T C^'fNj ( P\ - P 2 ) 


2uq(lu — ujq — | ir ) 


aL + FUt) 


(13) 


where u denotes the actual frequency of ajy. 

According to the definition, the real part of what inside the squaxe brackets is just equal to u). 
Thus we get an equation for u>: 


U — tUk + 


7TC 3 7 N i {P 1 -P 2 ) 
M (u> — LVq ) 2 + ^T 2 


(lx) — U>o). 


The expression of index of refraction follows immediately, 


u 

n = — = 1 + 
Uk 


7tc 3 7 

2 uJq (u — Uq ) 2 -f |r 2 


(u - UJo). 


(14) 


(15) 


The imaginary part of what inside the square bracket will be kj(u). The expression so obtained 
for Ki(u) is the same as eq.(7), except IV is replaced by T. The difference lies in that IV in eq.(7) 
already contains the contribution of Doppler broadening. 

We may generalize the above result by including the contribution of plasma free eletrons. 
Besides, since in the present case the wave number k is real , the operator aj^ decays with 
time, S +i kj must also do damped vibration accordingly. This means that u in eq.(14) must be 
analytically contiued to complex value which will be denoted by Q: 


— u>k + 


7TC 3 7 Nj(P !-P 2 ) uj 

2u% (Q - uq - |fr r ) 2(0. -if) 


(16) 


in which T is also replaced by Tr as in eq.(7). 

This is a third order equation for 0, it allows us to derive 0 = uj + iK for every given real value 
of Uk(uJk = kc ). We will choose the root whose real part is nearest to u k for photon. In the case 
of A = 5805 A, no problem appears. The so obtained curve uj-Uk is like the usual dispersion curve 
of gas. In this case the affection of medium polarization on 7 is indeed very small. 

In the case of A = 312A, the situation is different. For certain range of N e , the dispersion 
curve for photon and for collective atomic dipole (S+,kj) are totally mixed in the resonance region. 
This means the quantum of the polarization field has coupled to photon to form polaritons as in 
solidsi 6 ’ 7 !. In this case, the affection of medium polarization is significant. 

The occurence of such situation depends on the relative strength of coupling and damping. 
To show explicitly, let us examine the simpler case without contribution of free electrons. Now 0 
satisfies a second order equation. The two roots at u> k — uj 0 are 


0 


— + -V T ± 

4 



(17) 
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where 

a = - P *)> ( 18 ) 

just the product of the two couple constants between dj^- and S +ikj in eqs.(ll) and (12). 

When a < (weak coupling), the real part of the two roots given by eq.(17) equal each 
other, corresponding to the usual case as A = 5805 A. On the other hand when a > (strong 
coupling), the real parts of the two roots are different, corresponding to the situation of polariton 
formation. It is interesting to note that 


= r T 

2a «/( u>o) ’ 

therefore the two conditions < 1 and a » almost correspend to each other. Namely 

the region of polariton formation will cover the region of the pit discussed in last section. 
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Figure Caption 

1. Dependence of 7/7 on x for line 312A of ion Civ- x — W e /10 19 cm -3 , T = 5ev. 
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